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@ Single-Point Access to Data Distributed on Many Processors 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A description of the functions and data 
structures is defined that would be neces- 
sary to implement the Chapel concept of 
distributions, domains, allocation, access, 
and interfaces to the compiler for trans- 
formations from Chapel source to their 
run-time implementation for these con- 
cepts. A complete set of object-oriented 
operators is defined that enables one to 
access elements of a distributed array 
through regular aridnnetic index sets, giv- 


ing the programmer the illusion that all 
the elements are collocated on a single 
processor. This means that arbitrary re- 
gions of the arrays can be fragmented and 
distributed across multiple processors with 
a single point of access. This is important 
because it can significantly improve pro- 
grammer productivity by allowing the pro- 
grammers to concentrate on the high- 
level details of the algorithm without 
worrying about the efficiency and commu- 


nication details of the underlying repre- 
sentation. 

This work was done by Mark James of Caltech 
for NASA’s Jet Propulsion Laboratory. Further in- 
formation is contained in a TSP ( see page 1 ). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-42505. 


© Estimating Dust and Water Ice Content of the Martian Atmosphere 
From THEMIS Data 


NASA’s Jet Propulsion Laboratory, Pasadena, California 


Researchers at JPL and Arizona State 
University conducted a comparative 
study of three candidate algorithms for 
estimating components of the Martian 
atmosphere, using raw (uncalibrated) 
data collected by the Thermal Emission 
Imaging System (THEMIS) . THEMIS is 
an instrument onboard the Mars 
Odyssey spacecraft that acquires image 
data in five visible and nine infrared 
(IR) wavelength bands. The algorithms 
under study used data collected from 
eight of the nine IR bands to estimate 
the dust and water ice content of the 
atmosphere. Such an algorithm could 
be used in onboard data processing to 
trigger other algorithms that search for 
features of scientific interest and to re- 
duce the volume of data transmitted to 
Earth. 


The algorithms studied were based 
on regression models. In the study, 
the optical depths estimated by these 
algorithms were compared with opti- 
cal depths estimated in ground-based 
processing using fully calibrated data 
from both THEMIS and the Thermal 
Emission Spectrometer (TES). TES is 
an instrument onboard the Mars 
Global Surveyor spacecraft that also 
observes the planet at infrared wave- 
lengths, but at a lower spatial resolu- 
tion than THEMIS does. Of the algo- 
rithms studied, the one that 
performed best was based on a Gauss- 
ian Support Vector Machine regres- 
sion model. The test results indicated 
that this algorithm, operating on the 
raw data, had error rates that were 
within the uncertainty associated with 


the estimates obtained by the ground- 
based analysis of the fully calibrated 
data. This level of fidelity demon- 
strates that these algorithms are suffi- 
ciently accurate for use in an onboard 
setting. 

This work was done by Kiri Wagstaff Re- 
becca Castaho, and Steve Chien of Caltech 
for NASA’s Jet Propulsion Laboratory and 
Joshua Bandfield of the Arizona State Uni- 
versity. Further information is contained in 
a TSP (see page 1 ). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPQ-43590. 


© Computing a Stability Spectrum by Use of the HHT 

Unlike in the predecessor method, the mathematical sign of the damping is retained. 

Goddard Space Flight Center, Greenbelt, Maryland 


The Hilbert-Huang transform (HHT) 
is part of the mathematical basis of a 
method of calculating a stability spec- 
trum. This method can be regarded as 
an extended and improved version of a 
prior HHT-based method of calculating 


a damping spectrum. In the prior 
method, information on positive damp- 
ing (which leads to stability) and nega- 
tive damping (which leads to instability) 
becomes mixed into a single squared 
damping loss factor. Hence, there is no 


way to distinguish between stability and 
instability by examining a damping spec- 
trum. In contrast, in the present stability- 
spectrum method, information on the 
mathematical sign of the damping is re- 
tained, making it possible to identify re- 
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